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Abstract 



C^, We compute the distribution functions for gluons at very small x and 

not too large values of transverse momenta. We extend the McLerran- 
^ ' Venugopalan model by using renormalization group methods to integrate 

, out effects due to those gluons which generate an effective classical charge 

density for Weizsacker- Williams fields. We argue that this model can be ex- 
tended from the description of nuclei at small x to the description of hadrons 
. at yet smaller values of x. This generates a Lipatov like enhancement for the 

^ ' intrinsic gluon distribution function and a non-trivial transverse momentum 

dependence as well. We estimate the transverse momentum dependence for 
the distribution functions, and show how the issue of unitarity is resolved 
in lepton-nucleus interactions. 
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1 Introduction 

The problem of computing the distribution functions for gluons at very small x 
is an old one The gluon (and quark) distribution functions are computable 
in perturbation theory at large values of x, but at small x one encounters a so 
called Lipatov enhancement. The precise computation of this enhancement is 
subject to much uncertainty primarily because at some point in the evolution 
the density of gluons becomes so large that there are mutual interactions of the 
gluons which are ignored in the BFKL equation. In addition, the behavior at 
small X also involves knowing the distribution function at small Q^, and again 
non-pcrturbative information seems to be needed. 

Recently a different framework was advocated for the computation of the gluon 
distribution functions 1^ . The starting point of this approach is to view a hadron 
not as a collection of a small number of partons, but rather as a system with finite 
parton density. In the high density situation the natural way to describe the soft 
gluons is not as quasi- free particles, but as classical fields with large amplitude. 
These classical fields are generated by classical color charges which represent the 
valence partons. Once the high density effects are resummed into the classical 
fields one may apply weak coupling methods to calculate quantum corrections. 
In this approach, high gluon densities which prove so problematic in the BFKL 
context are a prerequisite for the description of the parton content of a nucleus 
wave function via classical gluon fields. 

A consistent separation in field and particle like degrees of freedom can be per- 
formed most easily in the infinite momentum frame. The object of computation 
is the intrinsic x and p_L contributions to the infinite momentum hadronic wave- 
function in the light cone gauge. Distribution functions are given as 

in terms of the intrinsic parton distributions as computed by taking the expecta- 
tion value of the number operator in the state of interest. 

In the infinite momentum frame the valence partons are strongly Lorentz con- 
tracted. If we then look for spread out gluon fields at a; ^ A^^^"^ interactions 
between valence partons and soft gluon fields eikonalize which indeed allows us to 
take the particle limit for the fast moving partons. Here A is the baryon number 
of a nucleus. For a hadron, a; <C 1. The hadron will indeed appear as an infinites- 
imally thin sheet on the scale of wavelengths associated with the momentum 
fraction x. (We will later see that we will have to regularize this source by giving 
it a large but finite momentum and a longitudinal extent of order i?/7 where R 
is its size in the rest frame and 7 is its Lorentz gamma factor. We will find that 
nothing in leading order of our computations depends upon the details of this 
regularization.) In addition, for a thick nucleus, since the number of sources of 
charge per unit area scales as A^^^, we may view the valence partons (quarks) as 
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classical color charges. Therefore, somewhat paradoxically the simplest problem 
to start with is computing the gluon distributions for a very large nucleus. 

We will find later that at very small x, the glue as well as valence quarks contribute 
to the charge density seen by a gluon. The giuons which contribute to the charge 
density are all giuons with an x larger than the x of the gluon whose structure 
function is being measured. Therefore the consideration discussed above for nuclei 
will apply to hadrons when at sufficiently small x so that the number of giuons at 
larger values of x is large. The advantage of nuclei is that large densities of charge 
are generated at larger values of x, and therefore lower energy per nucleon, than 
is the case for a single hadron. 

A solution of this problem would be useful in a variety of contexts. The approach 
we advocate involves knowledge of the nuclear wave function and is somewhat 
related to the approach of Mueller for heavy quarkonia |^ . Our approach in prin- 
ciple allows the resolution of various phcnomcnological problems which arise in 
the parton cascade model of particle production in heavy ion collisions Q . These 
models provide the initial condition for hydrodynamic calculations [||. A model 
which builds in the space-time structure we advocate and uses the information 
we have generated for the infinite momentum frame wave functions is given in 
Ref. §. 

The theory which results at the classical level is basically a Yang-Mills theory in 
the presence of a source 

J+ - S{x-)paix^) (1.2) 

The measure which generates the expectation values of gluon fields, corresponding 
to distribution functions is 



J [dA][dp] exp J (fx^^Trp"^ 



{x±_) expos' (1.3) 



where S is the ordinary gluon action in the presence of the external current J . 
The parameter /z^ is the valence color charge per unit area (scaled by a factor 
l/(iV^ — 1)). In leading order, the expectation value is given by a classical field 
which is a solution of the Yang-Mills equatiorj^ 

D^Ft'^^g'^r (1.4) 

which is then averaged over different values of p 

In the limit where the gluon field generated by these valence quarks is treated 
classically, the gluon field is a non-Abelian Weizsacker- Williams field, and has the 
form 

A+=A-=Q (1.5) 



^For conventions on the use of the coupUng constant please see the next section. 
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and 

A' ^0{x-)a'{x^) (1.6) 
The field a* 5 is a two dimensional "pure gauge" 

a' = -ic/V^[/t (17) 

i 

The physical justification for the non-Abelian Weizsacker- Williams field is that 
because the source of charge is confined to a thin sheet, the solution must solve 
the free equations of motion everywhere but on the sheet. The solution is therefore 
a gauge transform of zero field on either side of the sheet. The discontinuity of 
the fields across the sheet gives the charge density. 

It was suggested by McLerran and Vcnugopalan that this simple model should 
give a decent approximation for the soft glue distribution function. It turns out, 
however that the corrections to the distribution function calculated in this way are 
large at small x. Technically there are two sources for these corrections, although 
both have the same physical origin. 

First, as we will show the behavior of the correlation function calculated in this 
simple minded approach is singular at small p±. The fiaw in the treatment of 
Ref. 1^ was that the source of charge was not treated as an extended distribution 
which tends to a delta function only in the infinite momentum limit 1^, ^ . Phys- 
ically it is clear that the charge density is indeed spread out on the scale of the 
characteristic longitudinal momentum of the hard particles which generate this 
density. 

The second source of large corrections is basically the same small x enhancement 
as in standard perturbative calculations. As was shown in |^ the quantum cor- 
rections to the distribution functions calculated in terms of the classical fields 
become large at small x, the enhancement factor being the infamous alogl/x. 

The main goal of this paper is to show that the two corrections are physically 
related and to outline a solution to both problems. We will show that the small 
X enhancement arises from quantum fluctuations with large longitudinal mo- 
mentum. We show that such conflgurations may be successively integrated out 
by using renormalization group techniques reminiscent of the Wilson block spin 
method. This approach can also be interpreted in terms of the adiabatic or Born- 
Oppenheimer approximation extensively used in atomic physics. Integrating out 
hard quantum fluctuations is equivalent to including the harder gluons into what 



we call the charge density p in Eq. (L3) while calculating the distribution of 



the softer glue. It therefore leads to "renormalization" of the charge density and 
endows it with nontrivial and calculable longitudinal structure. 

This modified, momentum dependent distribution of source strengths leads to 
infrared non-singular correlation functions. We argue that the result is sensitive 
only to the average charge squared per unit rapidity per unit transverse area of 
the source. 
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The outline of this paper is as follows: 

In section ^ we study the classical problem of computing the fields associated 
with a source of charge which is extended in x~ . We find the general solution to 
this problem in light cone gauge. We compute the resulting distribution functions 
assuming that the source is random in x±_ and x~ , but with a weight of charge 
squared per unit rapidity per unit area which is specified. 

In section ^ we show by using the renormalization group techniques how to 
generate classical fields at some rapidity scale y. This involves perturbatively 
integrating out modes at larger values of rapidity (smaller values of x~). This 
integration generates an effective Lagrangian which has a self similar form, namely 
that at each step of the procedure it is similar to the McLerran-Venugopalan 
model but with a charge per unit area which is rapidity dependent. We show that 
this effective theory is equivalent to that discussed in section |^. We derive the 
renormalization group equations for the charge squared per unit area per unit 
rapidity as measured at some transverse momentum scale and rapidity y — 
yo + ln(xQ /a;^), where yo is the nucleus rapidity, and Xq R/l- The calculations 
in this section rely on several simplifying approximations, which we discuss. 

In section ^, we study the renormalization group equations for the charge squared 
per unit rapidity per unit area as a function of and y. This equation is closely 
related to the usual evolution equations for distribution function which appear 
in standard perturbative treatments. It can be viewed as a nonlinear version of 
DGLAP equation. We show that for p±^ much larger than the momentum scale 
associated with the charge squared per unit rapidity integrated over all rapidities 
larger than that at which we measure the structure functions (which we will refer 
to as xCj/jO^))! the nonlinearities in the RG equation become unimportant. In 
this regime the equation basically describes the double log DGLAP evolution. 
At lower momenta our equation can be thought of as a nonlinear variant of the 
BFKL equation, although to make the relation precise one would have to consider 
some virtual corrections in addition to those accounted for in our derivation. At 
low transverse momentum, we show that the evolution equation saturates. We 
discuss the consistency issues which are necessary for a solution of this equation 
within the set of approximations for which our derivation is valid. 

In the final sections we summarize our results. We show how our results are 
consistent with unitarity in deep inelastic scattering. We estimate the total cross 
section at fixed as x approaches zero. We argue that a computation of the 
charge squared per unit rapidity per unit area would allow a computation of the 
total multiplicity in hadronic interactions. We discuss the possible universality 
of our results and their possible generalization to the description of nuclcons at 
small X. We discuss some of the many problems which are not yet solved within 
the approach advocated in this paper. 
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2 Modification of the Source Strength Due to Ex- 
tended Structure in x~ 

In this paper we will use gauge potentials scaled such that the covariant derivative 
reads -D^[j4] = — iA^. The classical gluonic action is of the form — 
and hence a g^J^ term in the classical equations of motion. In this setup gauge 
transformations U are most economically parameterized via U{x) = exp [iA(a::)] 
transforming A as A ^ U [A — id] . We will be also using matrix notation, 
e.g. p = p"'T°', where T" are the normalized hermitian generators of the SU{Nc) 
group in the fundamental representation, 2TrTT^ = S""^. 

In the original McLerran-Venugopalan approach the source strength was as- 
sumed to have the form 

J+{x-,X^)^Six-)paix^_) (2.1) 

and to be distributed with the Gaussian weight 



[dp] exp 



d^xj_TTp\x^) 



(2.2) 



p^ 

where p,^ is the charge per unit area. 

The solutions to the Yang Mills equation in y4+ = gauge have vanishing A^ . 
Their transverse components A' are determined through 

V,d+A' + [A,,d+A'] = g^J+ (2.3) 

together with 

F'^ = (2.4) 
It was argued that the solution was of the form 

A'{x) =9{x-)a'{xi_) (2.5) 



In this solution, the commutator term in Eq. ( p.3[ ) was ignored since it involves 
the commutator of the field at the same point in x^ . 

Ignori ng t he commutator term is however not justified. It is clear that this term 
in Eq.( |2.3| ) is very singular and involves a product of S{x~) and 9{x^). To make 
sense of this structure, we must understand the evolution of the field across the 
delta function singularity. This can only be done if we know the structure of 
the source in x~ . In fact, as shown in ||] ignoring this problem leads to infrared 
singular distributions. 

Let us introduce the space-time rapidity variable 

2/ = yo + ln(a;o /x") (2.6) 

which will be useful for x~ > 0. We will assume that the source strength is 
non- vanishing only for positive x~ , and we will work in a gauge where the fields 
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vanish for x < 0. The rapidity j/o is the momentum space rapidity of the 
nucleus, and the parameter x^ is the typical Lorentz contracted size of the nucleus 
Xq - i?/7. 

In the next section, we will use renormalization group arguments to show that 
there is a non-trivial induced source strength extending beyond the volume occu- 
pied by valence partons which is driven by gluon modes at longitudinal momentum 
larger than that at which we measure the gluon distribution. This is parameter- 
ized by some strength of charge squared per unit area per unit space-time rapidity 
Q^), and by the charge per unit area at rapidities greater than y 

poo 

x{y,Q^) / dyy(y',Q^) (2.7) 
•'y 

The parameter appears because we must specify at what value of we 
are measuring the distribution function. It has precisely the same meaning as in 
perturbative QCD calculations, namely the transverse scale at which a parton is 
resolved . It should not be confused with the intrinsic transverse momenta of 
the fields. 

Accounting for the space-time rapidity dependence of the source strength, we 
therefore are lead to consider the distribution 



J [dp] exp 



(2.8) 



In this equation, p is the charge density per unit transverse area per unit space- 
time rapidity^. In the previous work we took great pain to argue that the charge 
could be treated classically on transverse scales which are large compared to the 
density of partons per unit area. This was because on this scale there is a large 
number of partons contributing to the source, and therefore the charges were in 
a large dimensional representation of the color group. This allowed a classical 
treatment. 

The longitudinal structure is a new ingredient. Why can we still approximate 
the partons (gluons) that couple to soft glue by a classical source? The physical 
reason is easy to understand: for these high momentum gluons the coupling is 
weak, so that to change the field (the soft glue) by a correction of order one, one 
must have many (hard) gluons contributing. In the next section, we will see that 
the induced source of gluons is in fact slowly varying in rapidity 



d^p dp 
dy^ dy 



2 "37 (2-9) 



Again, we will see this justified in more detail in the next section. 

^The parameter lJ?{y) controls the magnitude of the fluctuations of the charge density at 
fixed rapidity y. Since there is no charge density at rapidities greater than the rapidity of the 
r iucl eus j/Oi the function fi^(y) should vanish for y > yo. The rapidity integrals in Eqs. ( p.7[ ), 
!.8h are therefore effectively cutoff at this upper limit. 
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Therefore these sources of charge come from an extended region of space-time 
rapidity with a typical contribution at a rapidity far greater than that of the field 
we are computing. The source will therefore appear to be infinitesimally thin in 
the variable x~ . 

We must solve Eq. ( ^.3[ ) in the presence of source with a prescribed rapidity 
dependence. This is best done in terms of the rapidity variable y introduced in 
(PI). Eq. O) becomes 



Di—A'::^g^piy,xi_) 
dy 



(2.10) 



In this equation, because of the extended structure in rapidity, the term which 
involves the group cross product of two A'' fields cannot be ignored. 

The formal solution to this equation can be found by introducing the line ordered 
phase 



U{y,x±) = Uoo,y{x±) = Pexp 



I I dy'A{y',x_i_] 

V 



(2.11) 



representing a parallel transport operator along a straight line at fixed x~^ and 
x± connecting ?/ to cx) {x~ = 0). Recall, that due to vanishing of the transverse 
magnetic field {F''^ = 0), the vector potential should be a "two dimensional pure 
gauge" . We let therefore 



A\y,x^) = iUV'U-^ 
This leads to the equation for A 



(2.12) 



(2.13) 



The above equation may be solved directly numerically. Imagine we have a grid 
in rapidity y and transverse coordinates. We define the lattice spacing in rapidity 
as ay. The above equation can be written as 



V^A{y,x^) 



(2.14) 



= -.9^ P exp 



i I dy'A{y\x±) 

y+ay 



Piy,X-L) Pexp 



I I dy'A{y',x^) 

y+ay 



The solution at y depends only upon the function A at larger values of rapidity. 
This equation may therefore be solved iteratively starting at some maximum y^ax 
beyond which the source vanishes. 

It turns out, however that we do not need to know an explicit solution in order to 
calculate the distribution function. For that we have to perform the integration 
over the source strength p{y, x^_). Let us change the variables in the path integral 
Eq. (|2.§|) from p(y, x^) to A{y, x±^. The Jacobian of this transformation does not 
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depend on p. This is easily verified noting that the Jacobi matrix is triangular 
in y and due to the y orderings involved in the relation between p and A has no 
interactions on the diagonal. We find therefore that for any function 0{p) 



J [dp] exp ~ Jdy Jd\± 



a Trp(j/,x_L) 



0{P) 



(2.15) 



j [dK\ exp - j dy j rf^a;j_ 



Tr {V^K{y,x^)y 



0(A) 



Since the classical fields Ai are given as explicit functions of A, and our aim is to 
compute the distribution function 



Gtj{y,xj_;y',x'j_) = {A,{y,x^)Aj{y',x'j_)) 
this form of the path integral is very convenient. 
We first note that 



(2.16) 



OC 

= J dy' f/oo ,y'{x±) (W^A{y',x±)) Uy' ^oo{x±) 



(2.17) 



Now we perform the integrations over A by expanding the path ordered phases. 
It is most conveniently done by expanding the exponentials to first order on 
the rapidity grid with grid spacing Uy. This is a valid procedure as long as the 



function p^ 



is regular in the sense that lim ayp [y) — 0. We then perform all 



possible contractions with the propagator corresponding to the Gaussian weight 
in the path integral over A. Let us group together terms of the same order in the 
coupling constant. In zeroth order we have 



dy p\y,Q^)^iy'j=^ix±,x'^) 

V I V 



(2.18) 



Some comments arc in order concerning the inversion of the operator V*, since 
there is an infrared singularity in the inversion. Recall that the sources of interest 
ultimately arise from individual nucleons. Therefore all effects of sources die off 
at transverse size scales larger than 1 / Aqcd • The charge itself averaged over such 
transverse size scales also vanishes. This means that the Green's function should 
be defined with boundary conditions that ensure its vanishing at distance 1/Aqcd- 
In other words, whenever an infrared cutoff is needed for proper definition of an 
inverse of a differential operator, it should be taken of the order of 1/Aqcd- We 
will see that the quantities of physical interest are only very weakly dependent 
upon this non-perturbative length scale, but nevertheless such dependence does 
not disappear entirely. 



Here and in all that follows, we will define 



(2.19) 
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where 7(0) denotes an (infrared divergent) constant which ensures the vanishing 
of this Green's function as x± approaches the infrared cutoff 1 / Aqcd • Fortunately, 
the correlation function we will calculate below does not depend on the value of 
7(0) and the only infrared sensitivity that remains is through the logarithmic 
term in Eq. ( |2.19| ). 

In the first order, a quick computation gives 

1 r /"^ 

G\^{y,xr.y',x'^) = —Jab{-N,)g^ dy^^y" ,Q^) 

(V,V;.7(2:± - x'^)) [7(xi - xl) - 7(0)] (2.20) 

In this equation, is the number of colors. 
Similarly, in n'th order, we find 

(n+l) 



Gl.^,,{y,xr,y',x'^) = {-lT5ab 



(n + l)! 



d2/V(y",Q' 



V,V'7(x_L - x'^) {^{x^ - x'^) - 7(0)}' 



(2.21) 



The tadpole terms which take care of subtractions of 7(0) appear through the 
normal ordering of the path ordered exponential. This calculation can be found 
in Appendix A. 

We can now sum the series and find a representation for the correlation function 
as (assuming y > y') 



{y,xr,y\x'^) 



5'"' (v,v;7(2:± - xD) 



1 



N,[j{x^~x'^)^j{0)] 



(1 - exp {g^N.xiy, Q^Mx± ~ x'^) - 7(0)]}) (2.22) 



where x(y, Q^) is the total charge squared per unit area at rapidity greater than 
the rapidity y introduced in (pTTj). 

Finally, for the distribution function we get 

4(iV2- 



Ga 



Nax\ 



1 - (.xIA^qcd) 



(2.23) 



This correlation function has an amusing structure. At small transverse distances 
where 7 approaches zero, the correlation function tends to the perturbative cor- 
relation function, that is its value in lowest order in an expansion in . At short 
distances, the theory is perturbative. 

At large transverse distances (but of course still much smaller than 1/Aqcd), the 
correlation function dies off like It's Fourier transform at small momenta 

therefore behaves as 



G(2/,jy' = y,fci)^ln(fci/g47V^x(y)) 



(2.24) 
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Fig. 1: The distribution at fixed x as a function of intrinsic transverse momentum k\. 
We obtain considerable softening at small k\ compared to the perturbative l/k\ behavior 

This is in contrast to the behavior at larger transverse momentum where this 
correlation functions rises like 1/fc^ as fcj^ decreases, in agreement with the per- 
turbative result. The correlation function is therefore much softened at small k^. 
This behavior is shown in Fig.Q The characteristic momentum scale which differ- 
entiates between the non-perturbative and perturbative regions is k\ ^ g^Ncxiv): 
that is times the charge per unit area at rapidities greater than the rapidity 
at which the correlation function is measured^. This non-perturbative regime is 
nevertheless a weak coupling regime. Only when fc^ ~ ^qcd does the coupling 
become strong, and weak coupling methods can no longer be used. 

It is worth noting that the dependence upon Aqcd is very weak. At large trans- 
verse momenta the Fourier transform of the distribution function G does not 
depend on Aqcd- At large separations there is saturation, and there is again no 
dependence upon Aqcd ■ The dependence is really only in the region of very small 
momenta (fc_L <IC a^x): where our approximation is in any case not valid. 

This result is almost consistent with the structure which was advocated by McLer- 
ran and Venugopalan |^ . They had argued that at small transverse momentum, 
the above correlation function should approach a constant. It does up to loga- 
rithmic corrections. (The line of reasoning in |^ was however incorrect since it 
was based on an analysis of an equation that did not properly handle the induced 
charge associated with the gluon field, that is the [Ai, d^A^] term in the equation 

•^As we will see in the next section, the contribution of the gluons to x is proportional to Nc 
at large A^c. In the large Nc limit the coupling constant scales as g^Na = const. The crossover 
scale therefore has the correct large scaling behavior. 
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which determines the gluon field in terms of the external charge density). 



Renormalization Group Improved Charge Dis- 
tribution and Gluon Field 



In previous work, it was shown that radiative corrections to the distribution 
functions computed in the McLerran-Venugopalan model are large 0. The first 
order correction comes from the diagram, Fig.|^. 

The modification of the distribution function is of order ln(l/x) ln(fc^/as/i) 
at large k± and small x. For k± ~ a^/i, the natural k± scale in the problem, 
the corrections are of order asln(l/a;). In any case, for small values of x, these 
corrections become large and cannot be ignored. 



In this section, we will set up a renor- 
malization group procedure that sums 
up these corrections. The method of 
analysis is the following: We first con- 
sider the bare McLerran-Venugopalan 
model at a fixed valence charge den- 
sity, and a fixed ultraviolet cutoff in 
the longitudinal momentum . Phys- 
ically P"*" is of the order of the longi- 
tudinal momentum of the nucleus. It 
sets the scale of size for the longitudi- 
nal extent of the nucleus. 



A. 



A. 



I 



A 



A 



The renormalization group is imple- 
mented by considering the effective 
Lagrangian at a scale of momentum 
fc+ much less than P+ but where 
as\n{l/x) ^ 1. To generate this ef- 
fective Lagrangian, we integrate out 
quantum fluctuations with momen- 
tum fc"*" < g"*" < P"*". This procedure. 



Fig. 2: The leading perturbative correction 
to the gluon distribution beyond the clas- 
sical field approximation. The fat lines de- 
note gluons with large longttudinal momen- 
tum, p"*" . The momentum of the thin lines is 
fc"*" . The large logs come from the kinemati- 
cal region k+/p+ < 1. 

as will be seen generates a new effective Lagrangian of the same form as the origi- 
nal one, but with an additional charge squared per unit area. The typical scale of 
fluctuation of this additional charge squared per unit area is fj,^(jj,Q^)dy where 
dy = — ln(a;) and is a typical transverse momentum resolution scale at which 
the gluon distribution is ultimately measured. 

Since the form of the Lagrangian is unchanged under integration of these high 
momentum modes, except for the overall scale factor /i^, the procedure can then be 
repeated and yet lower momentum modes can be integrated out. Importantly, as 
long as the coupling constant is small and also ag ln(a;i/x2) is small the quantum 
fluctuations can be integrated out perturbatively, so that the computation is 
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controlled. This pcrturbativc treatment considers fluctuations around the classical 
solution in the region of the phase space which is being integrated out, as small. 
It is important to realize, however that the classical solution itself in this region of 
momenta is changed relative to the one at the previous step of the RG procedure. 
This is because it solves classical equations in the presence of the additional charge 
M^(2/) Q^)dy- This allows us to generate an effective Lagrangian at some scale of 
.T <C 1 . This is in a region where the naive McLerran-Venugopalan model would 
have broken down. 

Analogously, we can also allow the transverse momentum cutoff to be changed 
independently by a renormalization group transformation. This corresponds to 
perturbatively integrating out quantum modes in a phase space region with trans- 
verse momenta between Qi and Q2- This latter RG transformation is the counter- 
part in our approach of the standard perturbative renormalization group scaling. 

In the process of doing these transformations we develop a set of renormalization 
group equations for /^^(t/, Q^). These equations determine the rapidity and 
dependence of this parameter. 

We also determine the equation for the gluon field. It will turn out that this 

equation is a little more complicated than that in the previous section, since the 
induced charge depends upon the gluon field strength squared at the previous step 
of the renormalization group analysis. We argue that it should be a reasonable 
approximation to replace this fleld strength squared by its average value, in which 
case the equations described in the previous section can be derived. There are cor- 
rections to this approximation which are in principle computable. It is precisely 
this approximation which makes the fluctuations in the charge density uncorre- 
lated in space-time rapidity. Inclusion of these corrections will induce correlations. 
These correlations will however on the average not contribute to building up the 
charge density. 

It should be noted that there are other sources of correlation in longitudinal phase 
space. These arise from the classical field itself which is recomputed at each stage 
of the renormalization group analysis. Although it is true that the source of the 
color field is largely uncorrelated, for a given source, there arc still long range 
correlations built into the color field which would yield non-trivial multiplicity 
correlations in rapidity. 

Since the process closes under iteration, it is sufficient for us to show how we 
integrate the degrees of freedom as the Lagrangian changes scales in the Nth to 
the N-l-lst step of the renormalization grouping. This is what will be demonstrated 
below. 

We begin our analysis with the McLerran-Venugopalan action 




(3.1) 



where 



J+{x)=6{x-)p{x^) 



(3.2) 
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We are of course working in = gauge. 

Now suppose we have a solution to the classical equations of the form 

A+ ^ A- ^0 

A' = e{x-)a'{x±) (3.3) 

It is understood in the above expressions that the longitudinal delta function (as 
well as theta function) is regularized on the scale l/P^- Here is the typical 
longitudinal momentum of the fluctuations which have been integrated out in 
the previous step of the RG. At the very first step, of course P+ is the typical 
momentum associated with the nucleus or hadron. Clearly, the knowledge of the 
precise structure of the charge density on the scales of order 1 /P^ is necessary 
to determine the behavior of the classical solution at these scales. 

However, in the following we will only need to know the structure of the solution 
at longitudinal momenta much smaller than P^. This is so even though we will 
integrate over the fluctuations in the entire momentum range 

<k+<P+ (3.4) 

where P]^_^^i -C P^, but is still large enough so that as\n{P^ / P^^^) <C 1 (here 
as is evaluated at the scale x ^-nd x is assumed to be x ^ Aqcd)- The reason 
is that the dependence on the upper cutoff is only logarithmic, and the bulk of 
the contribution comes from much smaller momenta. To leading order therefore 
the results do not depend on the precise behavior of the classical field at the 
upper cutoff scale. Such is the magic of the logarithm, which enamored so many 
field theory practitioners. At momenta far below the cutoff the classical field does 



indeed have the structure A''(k) oc l/k'^a''{k±), which is equivalent to Eq. (3.3). 

To compute the effective action, we must integrate out the fluctuations around 
the classical solution. So long as we only generate an effective action at a scale 
fc+, so that as ln(P^^^/fc+) -C 1, then these fluctuations are small. We therefore 
have three types of fields to consider. There is the classical background field, the 
small fluctuation field at the scale of interest and the fields at lower momentum 
scale. We need only keep terms in the action which are at most quadratic in the 
small fluctuation field. We will denote these fields in the following manner: 

The field will be the classical background field. The label TV refers to the N'th 
step in the renormalization group procedure. This classical background field will 
be modified as the renormalization group procedure iterates. We will also write 

A^i = eNix')a'^ixA^) (3.5) 

In this equation, as mentioned before the step function is a step function on 
distance scale larger than that which we have previously integrated out. 

There are the small fluctuations fields at the step N which are within the mo- 
mentum range that we integrate out. We will refer to these fields as SA'^ 
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Finally there are the fields which are fluctuations around the classical solution at 
momentum scales much less than that where we perform the integration. These 
fields are not small. They are denoted as A'^ . 

The contribution to the effective action associated with the small fluctuation field 
is 

SS^^J d\^dx+^j {^-5A^{k+)D-^\k\p+)6A^{p+)^ + 

In this equation, the momenta fc+ and p"*" are in the range between the cutoffs 
Pn — l^^l: \P^\ — ^N-i- momentum is typically much softer than the 
lower cutoff, and therefore also much softer than The quantity is the 
inverse propagator in the background field. It depends on both the fields and 

We have approximated the linear term in the small fluctuation by keeping only the 
eikonal part of the interaction vertex, that is the coupling between the transverse 
components of the hard field and the minus component of the soft field. This will 
generate an effective action with only + components of currents affected by inte- 
grating out the high momentum modes. The terms we neglected are suppressed 
by factors fc+'/fc'*' and are sub- leading in the small x region. 



It is also understood that the transverse momenta of all the fields in Eq. (3.6) are 
bounded from above by some transverse cutoff Q. This is consistent with both, 
the BFKL approach, where all transverse momenta are roughly the same, and 
the leading log (or double log at small x) Altarelli - Parisi evolution, where the 
momenta are bounded by the momentum of the external probe. By imposing such 
a cutoff we restrict ourselves to transverse momenta which are not parametrically 
large. This point can be appreciated by examining the Feynman diagrams. Con- 
sider for example the diagram that gives the leading correction to the distribution 
function beyond the classical field approximation. It is depicted in Figj^. The cor- 
rections of this type with arbitrary number of insertions of the background field 
have been calculated in Q . For our present purposes it is enough to consider the 
classical field expanded to first order in the charge density p. The diagram then is 
precisely the same as that of the standard perturbation theory. After the integra- 
tion over the frequency k" is performed, the correction to distribution function 
is proportional to 



k+klJ Pi Jkt '^^ (p+ + fc+£l^ ^^'^^ 



The transverse momentum integration in this expression is cutoff not at the scale 
k\ (which is of the order Q^), but rather fc^/x, which at small a; is a very 
large scale. Physically the part of the integration region above k\ corresponds to 
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emission of jets which are much harder than the probe. Precisely the same problem 
is encountered in the standard perturbative treatment [jll] . These processes have 
to be considered separately, and at this point we will disregard them. 



Eq.(B.6) looks very suggestive. Introducing the notation 

Sp-ix)^2fabcafix)SAf^{x) (3.8) 

we see that the linear coupling term between the soft field and the hard fluctuation 
is of the form 

2TTdp{x)A-{x)d{x-) (3.9) 

We would therefore like to integrate in the path integral over those components 
of the fluctuation field SA, which are "orthogonal" to Sp. In other words we would 
like to change variables from 6Aib to Sp"" and some X"-, and integrate over X. 
In fact to get the result to the leading log accuracy it is not necessary to do 
it explicitly. Since 5p is linear in the fluctuation field, and the integral over the 
fluctuation is Gaussian, it is clear that the result of the procedure described above 
will be of the form 

[dSA] expjiJS'} (3.10) 

= M{a)J[dSp]e^pl^-J^ ^Sp^ix)Sp''iy)[Sx]Jix,y) + tSp-ix)A-^{x)^ 

Here M is the contribution of the determinant which arises in the Gaussian in- 
tegration over X. To the leading logarithmic accuracy this contribution can be 
ignored. This amounts to neglecting the loop corrections with all particles in the 
loop having the longitudinal momentum in the same range < p"*" < ^w-i- 
Corrections of this type do not give large contributions at small x |ll|| . It was 
also shown in the previous analysis that such contribution could be ignored 
at small x for modifications to the Weizsacker- Williams background field. 

The matrix Sxab{x, y) is given by 

5Xab{x,y) = 4z 1 ^r''f-ia'i{x^)a'^^{y^)Df^''f{p+,x^,x+,y^,y+) 

(3.11) 

To proceed further, we need to know the structure of the propagator of the hard 
fluctuations . Since the longitudinal momentum scale in the propagator is 
large we can use a no recoil, or eikonal approximation to incorporate the effect of 
interaction with the softer fields A'^~. The calculation is given in Appendix B. 

In the eikonal limit, the propagator is 

D^/\K+,z-;x+,x^,y+,yi_) (3.12) 
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S^jSix^ - y^f-^^^iO {signK+{x+ - y+)) 



xPexp 



In this equation, A^^^j and hence the path ordered phase is in the adjoint represen- 
tation. In the phase, the dependence upon x~ can be ignored since the function is 
slowly varying on scales of the typical size corresponding to l/fc+. Here we have 
taken as the momentum conjugate to the difference of coordinates x~ — y~ , 
and z~ = {x~ + y~)/2 ~ a;~ ~ a;~ is the average position associated with the 
field. 

The expression above accounts only 
for soft fields with longitudinal mo- 
Agi A(,| menta smaller than that of the fluc- 

; [ tuation 6 A. In terms of Fcynman dia- 

y., y_, grams this corresponds to summation 

U p of the diagrams of the type depicted 

on Fig.||. In fact, the propagator 
also depends on the background field 
'x1l{x±) and this dependence is im- 
portant in parts of the phase space. 
These contributions are of the type 
A A Fig.^. We will come back to this point 

and discuss the importance of these 
terms later. Temporarily, however we 
Fig. 3: Same as in Fig^ but with addi- will disregard them in order to make 

the discussion conceptually simpler. 




As for the soft insertions, the following 
remark is in order. The propagator de- 
pends only on the minus component of 
the vector potential. On the classical 



Same as in Fig.^ but with addi- 
tional insertions of the soft external field. 
These diagrams give corrections to the dis- 
tribution function upon contracting the soft 
legs. They therefore correspond to the virtual 
corrections and are higher order in fluctua- 
tion fields. 

solution discussed in the previous section, this component vanishes. It is therefore 
only the fluctuations of A~ around the new classical solution that contribute to 
D. Since these effects are higher order in the coupling constant, we will ignore 
them to this order. Again, these corrections too are important at low transverse 
momentum. This point will be addressed in section |^. 

With these approximations the fluctuation propagator becomes very simple. The 
fluctuation of the charge density is time (x^) independent and local in the 
transverse directions 



Sxabix±,y±) 



1 

2 



dvNS'ix^ - y^) f-^^f"-^ a'r{x^)atix^) 



(3.13) 



Note, that since our fields have a built in cutoff on the transverse momentum, 
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the a{x±)a{x±) actually should be understood as averaged on a transverse scale 
size (Px± ~ l/Q^- 

We now make the approximation 

1 



(3.14) 



A. 



A. 



i..v...C_J 



The averaging in Eq. (3.14) is over 
the distribution of p. We believe this 
approximation should be adequate to 
describe the RG flow of x, especially 
at large Q. The fluctuations of p are 
very short range in the transverse di- 
rection. On the other hand the fields a 
are slowly varying, its transverse cor- 
relation length being of order 
There is therefore very little correla- 
tion between p{x±) and a{x±) at the 
same point. Also, due to slow varia- 
tion of a{x±) in space, approximation 
of a'^{x±) by an x± - independent con- 
stant should be good with accuracy 
g^x/Q- Although this approximation 
is true on the average, it ignores some 
of the correlations which are built into 
the longitudinal structure. It would be 
very important to study corrections to this approximation or better yet to fully 
incorporate the structure of Eq. ( 3.13| ) in the solution to the problem. This is left 
for further study. 

We get therefore, that the change in the charge density is governed by the pa- 
rameter 

-dyN{a^)Q2 



A A 



Fig. 4: Same as in Fig!^ but with ad- 
ditional insertions of the hard background 
field. These diagrams are important in the 
region of momenta Q of order ftsX- 'S'ee the 
discussion in the next section. 



g^TT [N^ - 

_l 

g^TT {N^ - 



1) 



N 



dPk , 



(3.15) 



Here G°£(y, k±) is the transverse Fourier transform of the gauge field propagator 



Eq. (|2.23| ) . The variation of x due to the change of the transverse cutoff is also 
easily calculated 



N, 



1 



1 



N 



(N^ - 1) (27r)2 g^TT 



Y^dyp Gmy,Q) 



(3.16) 



The change in the effective action is therefore 

exp{i(5S'} — I [dSpN] exp ( — / d^x± 



-TTSp%{y,x±) 
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X exp i / d^x± / dx'^SpN{x±)Aj^{x±, x^) (3-17) 



Now we identify some typical space-time rapidity for our source with the momen- 
tum space rapidity. We expect that ^space-time ~ 2/mom- Let us define 

N 

y = yo-^dyi (3.18) 

i=l 

where the the right hand side is the momentum space rapidity shifts induced by 
integrating out the different scales. We will see that the left hand side has an 
interpretation of the space-time rapidity. The rapidity yo is the beam rapidity. 
We see that the equation for the evolution of x generated by our rcnormalization 
group procedure can be written as 

dx{y,Q^)/dydQ' = -J^-i-^ Gtny,Q) (3.19) 

This equation can be either formulated as a BFKL type equation, when one does 
the integral of Q first, or as a DGLAP type equation if one does the integral over 
y first |14[] . Notice that it is a non- linear generalization of both equations, since 
the right hand side of the equation is a function of x- 

This equation has a simple physical interpretation. The factor of Nc/{N^ — 1) is 
the charge squared per gluon. The number of gluons contained in our classical 
field is 

where S is the area of the nucleus. What our analysis has shown is that the change 
in charge squared is entirely due to the change in the number of gluons due to 
new phase space opening up. Of course this is a non-linear problem in general 
since the source of charge changes the classical background field in a non-trivial 
way. 



Let us also define 



and 



Piy,x±) = 5pN{xA_)/dy (3.21) 



A-^{x^,x+)=A-{x^,x+,y) (3.22) 

Now we must compute the change in the classical field. Since the change in the 
classical field is small, we see that if we write 



(3.23) 
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we can linearize the equations for SA'^. We find that SAj^^-^ d = ^ ^^"^ that 

D\A%)d+SA'i^, + [SA%,d+A%] = g^SN{x~)6pN/dy (3.24) 

In this equation, (5jv(a;~) means a delta function on the scale of our new classical 
field, that is as a regularized distribution it has its support on the scale 'i^/P^- 
Now we define SA'j:^ to vanish on the scale of the old classical field. Therefore only 
the first term survives. 

Upon identifying the index iV with the space time rapidity, we see that our clas- 
sical equation solves the equation we posited in the last section that 

D,iA'')^A^ =g^p{y,x^) (3.25) 
ay 

(The easiest way to sec this is to break up space-time rapidity into discrete in- 
tervals. Identify an index with each interval. Precisely the renormalization group 
equation for the field results.) 

We also see that the path integral measure for the fluctuating field is what we 
postulated in the previous section, with one caveat: We have omitted some contri- 
butions which on the average vanish in their contribution to the induced x- These 
terms generate non-trivial correlation in rapidity beyond which we compute. We 
will not further discuss their inclusion here except to note that they in principle 
are computable, and should be included at some point. 



4 RG Equations for xiu^Q' 



We now discuss the renormalization group Eq. ( 3.19| ). This equation determines 



how the color charge per unit area scales with rapidity and a transverse resolution 
scale size Q^. 

The consistency of our analysis requires that the solution to the renormalization 
group equation only involves information in the region where our approximate 
methods of computation are valid. It could easily happen that the region of inter- 
est in transverse space after several steps in the renormalization group procedure 
might drift to some value where our approximations are no longer valid. This 
might happen if at some rapidity y, the relevant typical values of became of 
order Aqqj-, or became much greater than y?{y,Q'^) where our classical source 
size approximation breaks down. It is plausible that the region of integration for 
the solution of the equations involves primarily the region of interest, since this 
is physically where the field originates, but we have no proof. In addition, the 
region of large where our classical methods no longer apply is probably cor- 
rectly treated even though the derivation above breaks down. In this region, the 
fields are weakly coupled, and our expression derived by classical means appears 
to be correct even in this region, to leading order in coupling. 
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The renormalization group equation may be formulated in the DGLAP form by 
first integrating over y as 



AT 11 fVo 

'^^Z^^' = (jV2 1 1) (2^)2 ^ I dy' G^ny\Q) (4.i) 

It may be written in the BFKL-Uke form as 



V 



dx/dy 



TTtAs^ r '^^'^ Gtt{y,Q') (4.2) 

1) (27r)^ g^TT Jo 



Note that this is a nonUnear equation. The dependence of its right hand side on 
X is determined by the solution of classical equations. 

At this point we want to return to discussion of the terms that we have neglected 



in deriving Eq. (3.19), namely the contributions of the diagrams of Fig.^, with 
insertions of hard background field. From the perturbative point of view those cor- 
respond to modifications of the gluon distribution due to mixing between the two 
particle and multi-particle (higher twist) operators. These diagrams can in prin- 
ciple be taken into account by using instead of the free propagator in Eq. ( |3.1l| ) 
the full propagator in the external field as calculated in |Q. 

Although this calculation has yet to be performed, it is easy to understand quali- 
tatively the main modifications it will bring about. First, even with the inclusion 
of the background field the additional charge density Sp will remain static. This is 
due to the fact that all the internal lines in the diagram Fig.^ have the frequency 
{p~ ) of the order of the on-shell frequency corresponding to the longitudinal mo- 
mentum < < P]^_i- It is much smaller than the on-shell frequency of the 
external line with momentum k~^. From the point of view of the emitted particle, 
therefore the coupling is always to the static source. The main effect of these 
extra insertions will be to modify the right hand side of Eq. ( 3.19| ) by adding 



to it terms nonlinear in (a)^. This effect however will be significant only for Q 
of order g^x- Physically, the diagrams of FigJ^ describe an emission of the soft 
particle with transverse momentum k± by a classical field ai{x). Clearly, as long 
as the transverse momentum of the emitted particle is larger than the inverse 
correlation length of the field, the particle is emitted locally. In this case the 
emission probability depends only on a'^{x) at the point of emission. In the local 
limit therefore the effect of these corrections will be of order a^(a;)/fc^. At large 
Q the main contribution to the distribution function comes from large fc^, and 
the correction due to nonlinearities is therefore negligible. At Q of order g^x ^^'^ 
smaller, the contribution of the diagrams in question is important. However, in 
the saturation regime Q ^ g^x they do not change the behavior qualitatively. In 
this region there is practically no running of x with y. The reason is, that since 
the correlation length of the classical field is of order {g^x)^^y the phase space for 
emission shrinks to zero at these values of momenta. 



The qualitative features of the solution of our RG equation are these. In the region 
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of large Q^, the equation approximately linearizes to become 



dy din tt 



X (4.3) 



This is precisely the double logarithmic approximation to the DGLAP equa- 
tion It would be solvable if it were not for the dependence of a on x- If 
we hold this fixed, we get approximately (assuming that x is a slowly varying 
function of y at some Qq) that 



X-^exp|^2^^y In QVg^j (4.4) 

The dependence of as upon x will of course modify the solution. 

In the region of transverse momenta in the vicinity of the crossover scale a^X 
the nonlinearities in the renormalization group equation become important. It is 
likely, that one of the effects of this will be that the transverse phase space in this 
region will be a very slowly varying function of Q^. It is then more convenient 



to turn to the form Eq. (4.2). Assuming this to be the main effect of the nonlin- 
earities, and approximating the transverse phase space by a constant P, we can 
write the solution as 

X = Xo exp {^^^Py^ (4.5) 

This has the BFKL - type behavior, growing as a power at small x. To calculate 
the value of the constant P we would have to include virtual corrections which 
have been neglected so far. 

Finally, in the saturation region where <^ ctlxiu^Q^)! the right hand side is 
constant up to logarithms. Here the solution is to a good approximation 

X - Xo + n{yo - y)Q^ (4.6) 

where k, is some slowly varying function. There is little change until (yo — y)Q'^ / Xo 
becomes of order 1. 



5 Unitarity, Total Multiplicity and Summary 

The issue of unitarity in deep inelastic scattering is related to the x dependence 
of 

f'^^ dN 

G{x,Q^)^ / dVy-T^ (5.1) 

Jq dxd^p± 

at fixed as x decreases. We have seen that at fixed p±, there are two separate 
regions for dN/dxd'^p±. The first is at large p^ ^ a^xi^)- In this region, the 
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integral above is xG{x,Q'^) ~ \n{Q'^)x{x) up to factors of logarithms of /i. As x 
decreases, this is a rapidly rising function of 1/x. 

At some point, for any fixed Q^, the parameter x will become < Q^. At this point, 
we are in the small p± region for the computation of dN/dx(Pp± . In this region, 
dN/dxcfp± ~ ln{p±) up to factors of ln(x) (It would be useful to determine these 
factors more accurately and actually compute the cross section in this region, but 
again this is beyond the scope of this paper.) Here the structure function xG 
has at most a logarithmic dependence upon x. There is therefore no obvious 
contradiction with unitarity. The dependence of xG on is also amusing, rising 
like some power of up to logarithms, until saturating at x- 

The total multiplicity produced in hadron-hadron collisions at x may also be esti- 
mated. Here we return to rapidity variables. On scaling grounds alone, the multi- 
plicity of produced gluons per unit area should be dNg/dynR^ xiv^V ~ !)• 
These gluons after production interact at high relative energy, and therefore 
largely elastically. The number of gluons should be approximately conserved. 
Later as a quark-gluon plasma is formed, the system expands approximately 
isentropically, so that the total number of gluons produced should be roughly 
the same as the number of pions. 

At this point, we do not have a full solution of the renormalization group equation 
in hand. Suffice it to say that one expects rapid growth of the parameter % as a; 
decreases at large p± . This should be much faster than a power of a logarithm of 
beam energy. The reason that this growth does not violate unitarity is because it 
is arising from an enhanced contribution at larger transverse momenta. 

The typical transverse momentum in this picture will go as the square of the 
multiplicity per unit area. The total deposited energy density at a typical forma- 
tion time t ^ 1/x will be of order x*. All of these functions are expected to be 
asymptotically somewhat rapidly rising functions of energy. 

To summarize: The results of this work are extremely suggestive. We have pre- 
sented a picture of low x gluon structure functions which has many of the intuitive 
features normally associated with the pomeron. The calculation presented here 
should be improved however in many aspects. 

We have made several drastic approximations in deriving the renormalization 
group equation. Let us once again point those out. 

First, we have neglected the virtual corrections. Those are generated by the dia- 
grams on Fig.^ when one contracts the external legs. Formally, as we mentioned 
in section ^, those are higher order in fluctuation and for that reason would seem 
to be sub-leading. However, some of these diagrams are known to contribute to 
the BFKL equation, and therefore must be important at least in some kinematic 
regime. This suggests that the generic form of the effective action which we have 



relied on, Eq. (3.17) is not quite complete. To see what is missing, let us con- 
sider for a moment fields in three ranges of the longitudinal momentum: the field 
A^(fc+), with k+ > the field B^(/+) with P+_-^ > 1+ > and the field 

Ci^i,{m'^) with m+ < P^. The integration over A and B generates the effective 
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action for C. Taking into account the soft insertions of Fig.||, means we should 
use for the fluctuations propagator the full eikonal expression Eq. ( |3.13 ) without 
setting the Wilson line factor equal to unity. The integration over the fluctuations 
of the field will generate the effective Lagrangian 



i / c?^xj_Tr5pjv_i(a;_L)- 



XN-i{B-+C-) 



(5.2) 



where 



and 



XN-i{B- + C-) - SxN-iW{B- + C-) 



W{A-) = Pexp 



(5.3) 



(5.4) 



In the next step the integration over the fluctuations of B leads to the effective 
Lagrangian for C 



(fx± / dx^ [dpN-iix±) + SpNix±_)]C {x±,x~^) 



I 



+i j d^Xi_Ti5pM{x±_) _ j„_. SpN{x±) 



6pN-l{xj_) 



(5.5) 



In the second term the brackets denote averaging over the fluctuation of the field 
B. 

iXN^B- + C-))b = ^XA^-i ![d5B] W{B + C) - (1 - -i)5xN-iW{C) (5.6) 



We see therefore, that the integration over the fluctuations in the range of mo- 
menta between P^_i and -P^, not only generates the additional charge density 
Spff, but also modifies the fluctuation amplitude of the charge density p which is 
generated by the higher momentum modes fc+ > P^_i- This modification is only 
important for the coupling of the density p to the fields with momenta m"*" < , 
since only in this case the longitudinal phase space is large and the correction fac- 
tor 7 is proportional to log(l/a;). In terms of Feynman diagrams this calculation 
corresponds to the virtual corrections of Fig. 5| and the factor 7 is directly related 
to the so-called non Sudakov form factor Jll|. 

It is important to note, that even though the form of effective Lagrangian Eq. ( |3.5[ ) 
is not precisely the same as considered in section ^, on the classical solutions where 
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A~ — the two indeed coincide. The solution considered in section || is therefore 
still applicable to the modified Lagrangian. The net eS'ect of the virtual corrections 
is to modify the running of the effective charge density xin) through the change in 



the right hand side of the renormalization group equation Eq. (3.19). This effect 
is calculable, and should indeed be calculated, but this is beyond the scope of this 
paper. We will only note, that since these virtual diagrams do not play any role 
in the perturbative double log DGLAP treatment, and our renormalization group 
equation reduces to it in the limit of large Q, we expect these extra corrections 
to be important only in the nonlinear regime Q ~ ctsX- 

As a second approximation we ne- 
glected the insertions of the hard back- 
A ground field. This was discussed in the 

previous section, where we have ar- 
gued that these corrections are also 
unimportant at large Q. Again, in 
principle these corrections are calcu- 
lable if we employ the full fluctuation 
propagator in the background field, as 
calculated in 0. 

The third approximation in arriving 



T 




1 
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o 

B 













Fig. 5: Virtual corrections of Eq. (5.1). 
The thickest line denotes the propagator of 
the field component A with the largest longi- 
tudinal momentum. The thick line that com- 
pletes the loop is the propagator of the field 
B. Finally, the thm external denote the ex- 
ternal field C , for which the effective La- 



ta Eq. (3.19) was to replace the square 
of the classical field by its average. 
This lead to the absence of correla- 
tions in rapidity for the density fluc- 
tuations. This approximation is also 
expected to be good for large trans- 
verse momenta, where the emission of 
soft fleld is local in transverse coor- 
dinates and therefore the additional 
density Sp is practically uncorrelated 



grangian is being computed. 
with the density coming from higher rapidities. At transverse momenta of order 
of the inverse correlation function of the classical field this approximation should 
break down. 

Clearly, the treatment of the nonlinear region Q < a^x in the present paper is very 
rudimentary. It is in great need of improvement, and we intend to address this 
problem in future work. In fact, at small transverse momenta the very notion of 
the momentum cutoff Q is highly questionable. The density - density correlation 
which is given by the diagram Fig.^, in momentum space is a very slowly varying 
function of the transverse momentum at k± ^ c^sX- K depends on the momentum 
logarithmically. In this range of momenta one can therefore approximate it by a 
constant. In our calculation precisely this is achieved by introducing the cutoff 
Q and using the eikonal approximation for the propagator D^, which results in 
a local correlator of density in the transverse coordinates. At momenta of order 
ttsX, however the correlator changes rapidly. Approximating it by a constant 
with some transverse cutoff should therefore result in an error of order one. A 
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more careful treatment will bring about nontrivial transverse correlations of the 
charge density. One should therefore expect, that an improved treatment of low 
transverse momenta region will modify the distribution for density fluctuations 
such that nontrivial transverse as well as longitudinal correlations will appear. 
The effective Lagrangian which would generate the classical solutions will be of 
the form 

J dydy'(fxj_(fx'j^TTp{y,x±)[fi'^{y,y',xj^,x'j^)] ^ p{y' ,x'^) (5.7) 

In fact, in a general case there is no reason to expect that the weight will be 
Gaussian, so that the weight function p, could itself depend on p. 

It remains to be seen how large in fact will be the effect of these improvements. 
We believe, that although quantitatively it should be significant, the qualitative 
picture presented here will be confirmed. If the above issues can be answered in 
a satisfactory way, then one can proceed to a detailed computation of hadronic 
interactions within this approach. Deep inelastic scattering from nuclei and Drell- 
Yan production in heavy ion collisions might be computed. The initial conditions 
for nucleus- nucleus collisions might be found in detail. The fluctuation spectra 
and correlations between fluctuations would also do much to verify the above 
picture. It would be very useful to have data on the structure functions directly 
by using nuclei in HERA. 
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A Normal Ordering the Distribution Function 

We have to calculate 

G,,(y,x^;y',xl) = J [dp] e^p J dy'^d^x^^-^^^j^^p^y'' , x'^)^ 

X t^U{y,x^)^,UHy,xi_)Uiy\x'^)V,U\y',x'^) (A.l) 
This can be cast in the form 

G,j{y,x±;y',x'j^) (A.2) 
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= ijdy' {Uoo,y' (V^A(y')) f/y.oo) {xi_) J dy' {U^^y, (V'A(y')) f/^',-) (xi))A 
y y 
with the correlation function 

(A(y,Xi)A(y',.xl))A = g^^\y,Q')S{y-y'hix^-x'^) (A.3) 



The function 7 is given in Eq. (2.19) and as discussed in section g is infrared 
singular. The leading dependence on the infrared cutoff resides in the constant 
term 7(0). Fortunately, all terms containing 7(0) cancel in the expression for the 
correlation function. To understand how these cancelations work in (A.3) let us 
first consider the normal ordering of individual link operators first. To do so, let 
us break up any link operator from ?/ to 00 into "infinitesimal factors" 



lim T\UoO,y,{XjL)Uy^^y^_,iXl_)- ...■Uy,,y{Xl_) (A.4) 



In the large k limit each Uy^^y^_-^ covers an infinitesimal piece of the total path 
with a fixed length A = i/m — ym-i- Due to the locality of ( A.3 ) in y it is clear 
that there will be no contractions between different factors in this product. For 
an individual factor however we may expand and perform the normal ordering 



ym yr. 



1 + i JdyA{y,x±)+i^ JdyJ 



3\ 



I- I dy J dy' A{y,x±)A{y',xj_) 

ym-i y 



+0(A^) 



1 + i dy A{y,xi_) : + 



-O(A^) 



J/m-l 

2\ 



ym-i 



(A.5) 



where we have kept all terms up to order A. This is the only non-suppressed 
tadpole contribution if the function fi^{y) is finite. As a consequence, we have 

Uy^,y^-i{x±)Uy^^,,y^_^ix±) 

ym / 4 \ 

= + * \dy A(y,a:^):+1 (^_ g^^^W j ^ dy ^,\yM'') 



+0(A2 



ym,ym-2 



(x_l) : exp 



dyi?{y,Q'') 



Vm-I 



(A.6) 



which immediately carries over to Uoo,y{x_\_) upon insertion into (A.4). The dan- 
gerous tadpole contributions therefore can be factored out from a link operator 
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by writing it in the normal ordered form. Using this result we find 

Gij{y,x±;y',x'j_) 



CO oo 



(A.7) 



= Jdy' J dy' (: {Uo.,y' {V'Aiy')) Uy,,^) {x±) : 
y V 

X : [U^^y, (V^A(y')) C^s'.oo) :)a 



X exp 



oo oc\ 



OC 'DC 

jdy' jdy' 

V V 



Gij{y',x±;y',x±) 



X exp 



+ J I dy' f,'{y',Q^) 

-y' y' 



oo oo\ 

J + J\dy' n\y',Q' 

\y' y' I 



The expectation value of the product of normal ordered fields we dubbed G 
does not contain any contractions within the individual /7's. This is now easily 
evaluated order by order and then resummed. Expanding the C/'s to zeroth order, 
we have 

Glf{y\ xr, y', x^) = S"'6iy' - y') g^p^iy' , Q^)ViVj^{xi. - xi_) (A.8) 

In the first order, a quick computation gives 

G'^\y',xr.y',XA_) (A.9) 

= G^^{y',x^;y',x^) {-g^){-N,h{x^ - x^) Jdy^^y" ,Q^) 

y' 

In this equation, Nc is the number of colors. 
Similarly, in n'th order, we find 



Gtt''{y',xr,y',xi.) = 



{g^YN^-f{x i. - XA 
n! 



dyyHy",Q') 



by' 



(A.10) 



X Gt^°{y',x^;y',x^) 

Summing up, multiplying the tadpole factors and then performing the remaining 
y' and y' integrals now allows us to write an explicit expression in which the 
leading infrared divergence cancel. Assuming y > y the result is 

G'^^{y,xr,y',x'^) = -S->'{v,V^jix^-x'^))-— ^— — 

(1 -exp{9^AfeX(2/,Q')[7(a;± -^Y) -7(0)]}) (A.ll) 
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B The eikonalized propagator 



where we have defined 



/•oo 

x{y,Q^)= / rft/V(j/',Q2) (A.12) 
Jy 

The quantity x(j/, Q^) is the total charge squared per unit area at rapidity greater 
than the rapidity y. 



B The eikonalized propagator 

In this appendix, we will derive an expression for the vector field propagator in 
the eikonal approximation. We will solve the equation of motion for the hard 
fluctuation field in the presence of external soft vector potential. Throughout this 
analysis we neglect the effects of the classical background field. 

We start with the transverse component of the Yang-Mills equations: 



D^F^"' = (B.l) 



and write the total field as 



where SA^ is the hard field describing the fluctuations with high longitudinal 
momentum and is the soft field with small longitudinal momenta only. 

We assume that the only large momentum in the problem is the longitudinal 
momentum of the hard field SAi. Therefore, in the equation of motion for the 
hard field we keep only those terms which involve derivatives of the hard field with 
respect to x~, the coordinate conjugate to large momentum p+. The equation of 
motion for the hard field then becomes: 

d-d-5Ai - i[s-,d-SAi] = D-{s)d-6Ai = (B.2) 

To calculate the propagator we need to find eigenfunctions of the operator 
D~{s)d-. In order to do this we write 

6A^{x) = e'''^SA^{x) 

where the eigenvalue X = and SAi is a slowly varying function of a;~. Then 
eigenvalue equation becomes 

D-{s) 6Ai{x) = 

which has the solution 



dz~^s~ {z'^ ,Xt) 

-oo 



(B.3) 
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where a is the color label, A is the eigenvalue index, el^-* and U(^a) a,rc the po- 
larization vector and color basis vector respectively. The eigenfunctions therefore 
are 



Pexp(— i / dz~^s {z~^^x ,Xt) 

J — oo 



The propagator is constructed as 



G'}';{x,y) 



dX 



5At'\x)6Af'\y) 



X — ie 



Since the soft field s is a slowly varying function of x~ , we can neglect its variation 

with x~ and write it as a function of z~ where z~ ^ {x~ + y~) is the average 
(— ) coordinate associated with the soft field. The expression for the propagator 
can be written as 

f dn+dn~ p-'i-P'^ -v~) p-ip~ (x* -y^) 



Pexp(— i / dz'^s {z~^,x ,Xt) 



(B.5) 



- ab 



Here we have used 



and 



The integration over p~ is straightforward and gives a factor proportional to 
9{x~^ —y~^) for positive p+. To get the propagator in momentum space, we Fourier 
transform with respect to relative and center of mass coordinates, x~ — y~ and 

z^ ~ {x^ +y~) to get: 



G1^{K+,k+,x+,y+,Xuyt) = lSij6\xt-yt)9{x+-y+)^ 



Pexp(— i / dz~^s {z'^,k^,Xt) 



(B.6) 



- ab 



where and fc+ are the momenta conjugate to {x~ — y~) and {x^ + 2/ ) 
respectively. For negative we get the same expression above with the argument 
of the theta function switched around and a relative minus sign. 
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